Surface Specific Heat of 3 He and Andreev Bound States 
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High resolution measurements of the specific heat of liquid 3 He in the presence of a silver surface 
have been performed at temperatures near the superfluid transition in the pressure range of 1 to 29 
bar. The surface contribution to the heat capacity is identified with Andreev bound states of 3 He 
quasiparticles that have a range of half a coherence length. 

PACS numbers: 67.57.-z, 67.57.Bc, 74.45.+C 



Unconventional pairing superfluids and superconduc- 
tors are sensitive to quasiparticle scattering at sur- 
faces since all forms of scattering are inherently pair 
breaking [[J. Depending on the boundary conditions, 
whether scattering is specular or diffuse, and depend- 
ing on the specific quantum state, the order parame- 
ter can be significantly suppressed. Correspondingly, 
quasiparticle bound states extend from a surface a dis- 
tance approximately equal to the coherence length of 
the bulk superfluid. These states were first discussed 
by Andreev[2j in order to understand the difference be- 
tween charge and thermal transport at superconducting 
interfaces, and they have been extensively investigated in 
unconventional superconductors. For example, the zero- 
bias conductance anomaly in tunneling experiments 
has been ascribed to low-energy, surface bound states and 
provides a key indicator of unconventional pairing0,[3- 
Andreev scattering and Andreev bound states (ABS) 
are essential characteristics of thin superfluid films pj of 
3 He and they dominate the properties of superfluid 3 He 
contained in the porous medium of silica aerogel[3 0- 
In the latter case, the ABS lead to gapless superfluid- 
ity as has been determined from their influence on heat 
capacity [IJj and thermal conductance |Tl| . The bound 
states affect physical measurements that use probes such 
as vibrating wires0, for viscosity and thermometry 
experiments, and crystal oscillators [l^, [lf| for the mea- 
surement of the acoustic impedance. 

Recently Vorontsov and SaulsQ have calculated the 
contribution to the free energy and specific heat of An- 
dreev bound states in thin films of 3 Hc in the A-phasc. 
For a film on a solid surface having diffuse scattering 
boundary conditions they find that there is a suppression 
of the superfluid transition, T c , as well as a substantial re- 
duction of the heat capacity in the superfluid state near 
T c . Even for thick films, where the suppression effect 
on the transition temperature is negligible, the reduction 
of the heat capacity near the transition can be remark- 
ably large. In this letter, we present measurements of the 
contribution of these bound states to the heat capacity of 
superfluid 3 He near a silver surface close to the transition 
temperature. 

Previous experimental work on 3 He in confined ge- 



ometries, has taken one of two approaches. The first 
is to investigate 3 He thin layers, for example, films 
with a free surface for studies of superfluid density 
flowfr^. third sound[l^|. or in slabs having confinement 
on two sides, as was the case for a number of NMR 
experiments [l9(. The other method is to determine 
the effects of surfaces on 3 He constrained in a porous 
medium, with the corresponding advantage of a larger 
effective surface area. If the pore structure is larger than 
the superfluid coherence length, the system can be ap- 
proximated as a collection of randomly oriented planar 
surfaces. For measurements of the heat capacity this lat- 
ter approach is preferable. Earlier experiments |2fl \2 11 122| 
of this kind show that the heat capacity differs from that 
of the bulk, without a consensus for interpretation. Grey- 
wall suggested that there is a healing length of the super- 
fluid at the surface |20|. Others have argued |2fl that 
there is a broad distribution of transition temperatures 
of disconnected superfluid regions. In our experiment we 
use a high resolution temperature sweep method that can 
provide sufficient detail to explore the temperature and 
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FIG. 1: Heat capacity of both bulk and confined He 
obtained from a slow warm up trace at 11.31 bar. The 
solid trace is the heat capacity expected for bulk 3 He 
determined from Greywall's measurements]^]. The data 
points are our measurements. 
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pressure dependence of the heat capacity, and we keep 
some bulk helium present in the calorimeter as a ref- 
erence. Near the transition we observe a deficit in heat 
capacity with respect to the pure superfluid as was found 
in the thin film calculations of Vorontsov and Sauls. The 
model we develop is based on a surface specific heat from 
surface Andreev bound states. We find that the model 
can consistently account for our results as well as those 
from the earlier work. 

Our measurements were performed with the calorime- 
ter described by Choi et al. 01 f° r high resolution mea- 
surement of the specific heat of superfluid 3 He in sil- 
ica aerogel. There are three regions of 3 He inside the 
calorimeter. The first is the interior of the silver heat ex- 
changer constructed of sintered silver particles of micron 
size and has a volume V\ — 0.56 ± 0.01 cm 3 and surface 
area 2.9 ±0.1 m 2 . The second region is the open volume 
for bulk helium, V 2 = 0.29 ± 0.04 cm 3 . Finally, from our 
earlier studies|10j. we have a disk of silica aerogel with 
pore volume of V3 = 1.06 ± 0.01 cm . The 3 He in the 
volume V3 remains in the normal Fermi liquid state for 
all of the experiments reported here. We have previously 
determined its volume and heat capacity to an accuracy 
of 2%. We have subtracted this contribution, plus the 
calorimeter background, from our measurements and do 
not discuss them further. 

The samples were cooled by adiabatic demagnetization 
of PrNis and the calorimeter was isolated from this re- 
frigerator with a superconducting cadmium heat switch. 
The temperature of the sample cell was measured ev- 
ery thirty seconds using a SQUID based mutual induc- 
tance bridge for measurement of the magnetic suscep- 
tibility of a paramagnetic salt, La diluted CMN. Once 
the cadmium superconducting heat switch was open, the 
sample cell warmed at a rate, T, from an ambient heat 
leak, Q, typically 0.1 nW. Occasionally we applied exter- 
nal heat pulses to check consistency and to calibrate this 
heat leak. Then the heat capacity was determined as, 

dQ = dQdt_ = Q_ 

~ dT ~ dt dT~ f 

The advantage of using slow warming traces over the adi- 
abatic heat pulse method is higher resolution. A heat 
pulse typically causes a temperature jump of 50 ~ 100 
fiK. In a slow warm-up trace, the temperature change 
for each point is less than 1 /xK. However, such a small 
signal inherently results in poor signal-to-noise in deter- 
mining T . This can be overcome by averaging adjacent 
data points provided that the warm-up rate is adequately 
slow and stable. We used averaging to smooth the data, 
thereby decreasing our temperature resolution to 10 /iK. 
All of our slow- warming data is reproduced by our pulsed 
heat capacity measurements, albeit with lower resolution 
in temperature . 

On cooling through T c we observe a sharp, resolution- 
limited, increase in the heat capacity, shown in Fig.l, 
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FIG. 2: 8C — (C — C s ) is the difference between the mea- 
sured heat capacity, C, and that of the bulk superfluid, 
C s , as a function of temperature at 11.31 bar. The inset 
is a sketch of the volume distribution in the calorimeter. 
Vi is the fluid inside the silver heat exchanger and V2 the 
volume outside. In our model, Andreev bound states re- 
side within a distance a£(T, P) from the surface in the 
volume Vi; the rest of the helium in Vi and all of that in 
V2 is taken to be bulk superfluid. The model calculation, 
given by the smooth curve with a constant scale factor 
q = 0.48 ± 0.08, agrees well with the data. 

followed by a smooth increase and then a decrease over a 
range of temperature. For reference we directly compare 
our results in this figure with the heat capacity measure- 
ments of bulk superfluid 3 He performed by Grevwall|2fj|. 
The central question we address is, what is the origin of 
the difference between these results. For bulk helium we 
know that the jump in heat capacity at T c , AC S (T C ), cor- 
responds to that of a BCS pairing system, enhanced by 
strong coupling 23] . The heat capacity then falls rapidly, 
approximately proportional to T 3 . Consistently, in our 
data we find that at T c there is a sharp increase in the 
heat capacity and it is natural to identify this jump with 
the bulk helium in our calorimeter. In Fig. 2 we show the 
difference between the measured heat capacity and that 
of the bulk superfluid for the same volume, SC — C — C s 
as a function of temperature. The magnitude of the dis- 
continuity in SC at T c corresponds to the amount of he- 
lium in the silver heat exchanger given by the volume 
ratio Vl/(Vi + V 2 ) = 5C{T C )/ AC S {T C ) and is plotted in 
Fig. 3. The apparent volume Vi, deduced in this way, 
is 0.40 ± 0.02 cm 3 . As expected, it does not vary with 
pressure. The magnitude of the apparent volume is qual- 
itatively consistent with an independent measurement, 
Vi = 0.56 ± 0.01 cm 3 . Apart from experimental uncer- 
tainty this discrepancy reflects difficulty in making an 
accurate extrapolation to T c , which we discuss in greater 
detail below. 

Below T c the behavior of the heat capacity must be 
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attributed to the combination of the surface dominated 
heat capacity in the silver heat exchanger in addition to 
that of the bulk. The formation of surface bound states 
corresponds to transfer of spectral weight from above the 
energy gap to low energy (near the Fermi energy) with a 
density of states of low energy excitations of order that 
in normal helium. Their spatial extent from the surface 
is expected 0, 0| to be approximately that of the co- 
herence length given by £(T,P) = £o(P)(l - T/T^ 1 / 2 
where £o(-P) = tuup /27rfcsT c . Here Vp is the Fermi ve- 
locity and ks is the Boltzman constant. On this basis 
we propose a simple model where we take the surface 
heat capacity to be proportional to that of the normal 
fluid but constrained to a volume that scales as A£(T, P) 
where A is the area of the silver surface. Consequently, 
we write the surface contribution to the heat capacity as 
a(T,P)£(T 7 P)Ac n and we investigate the temperature 
and pressure dependence of the scale factor a. A picto- 
rial representation of this model is sketched in the inset 
of Fig. 2. The corresponding heat capacity is, 



C = a(T, P)£(T, P)Ac n + {Vi + V 2 
= (Vi + V 2 )c s + a(T, P)Z{T)A{c ri 



a(T,P)£(T,P)A)c s 
- c.) (2) 



and SC = a{T,P)i(T,P)A(c n - c s ). The scale factor, 
a(T,P) is the only unknown parameter necessary to de- 
scribe the surface heat capacity and expresses the tem- 
perature dependence and pressure dependence of the sur- 
face heat capacity beyond that given by the coherence 
length and the normal specific heat. 

We have used this model to interpret our measure- 
ments for various temperatures and pressures. For any 
given pressure, we find a is constant over the available 
range of temperature, down to T/T c w 0.7. In Fig. 2. 
our measurements of 8C at a pressure of 11.31 bar are 
compared with a fit to Eq.2 taking a to be temperature 
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FIG. 3: Measurements of the heat capacity discontinuity 
at T c interpreted as the volume in the heat exchanger, Vi 
for various pressures. 
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FIG. 4: The scale factor a for the surface heat capacity 
as a function of pressure. The present measurements are 
shown as open circles and are pressure independent with 
an average of a — 0.51. There is good agreement with 
earlier work from Greywall 20] (diamond) and Kishishita 
et al. |22jl (solid circle) interpreted in terms of our model. 



independent. The good agreement between the data and 
the calculated curve confirms that a is a constant with 
a best fit value of 0.48 ± 0.08 at this pressure. The sig- 
nificant down turn in 5C near T c in Fig. 2 is due to the 
strong temperature dependence of the coherence length. 
We have made this comparison at all pressures and the 
results are presented together in Fig. 4. The scale factor- 
appears to be both temperature and pressure indepen- 
dent with the average value a = 0.51 ± 0.15. In the con- 
text of our model this means that the spatial range for 
surface excitations, that we associate with quasiparticle 
bound states, is a half of a coherence length. 

Greywall 20] allowed for a healing length of superfluid 
3 He near the silver heat exchanger surface in his mea- 
surement of specific heat and he assumed its tempera- 
ture dependence to have the form (1 - (T/T c ) 4 ) -1 / 2 . We 
have reanalyzed his data with our model as well as the 
work of Kishishita et al.^j^. Both results are plotted 
in Fig. 4 where they are compared directly to ours. It is 
noteworthy that the silver sinter used by Kishishita et 
al. had an area-to- volume ratio of 12 x 10 6 m -1 , the 
one in the Greywall experiment was 3.4 x 10 6 m , and 
these can be compared with ours, 5.2 x 10 6 m _1 . There 
is excellent agreement among the experiments, although 
they are performed in a range of pore structures with 
area-to- volume ratios spaning a factor of three. This im- 
plies that different structures among the silver sinters do 
not play a role. The overall consistency of the data with 
the model, including the variables of pressure, tempera- 
ture, and different silver surface structures, provides com- 
pelling evidence that we are measuring a surface contri- 
bution to the heat capacity, rather than the heat capacity 
of disconnected regions of superfluid with a distribution 
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of transition temperatures. However, the model will not 
be correct close to the transition temperature where the 
coherence length diverges with increasing temperature 
approaching T c . There is a point, nominally a few percent 
lower than T c , where a£(T,P) = a£ (P)(l - T/Tc)" 1 / 2 
reaches V\/A = 193 nm. In our model, a£(T, P)A rep- 
resents the volume of the surface bound states and, at 
this temperature, they would fill the silver exchanger of 
volume V\. At 1.024 bars this excluded temperature re- 
gion is 3% of T c ; it decreases with decreasing coherence 
length at higher pressure. Additionally, there is a small 
suppression of the transition temperature for helium in 
restricted geometry. These effects may account for diffi- 
culty in extending our model close to T c and the corre- 
sponding systematic error from such an extrapolation in 
determining V\ as is shown in Fig. 3. 

From a theoretical perspective, quasiparticle scatter- 
ing at the surface is responsible for a non-zero density of 
states at the Fermi level which should give a heat capac- 
ity that is linear in the temperature in the low temper- 
ature limit. Our model for the surface specific heat has 
this temperature dependence at low temperatures where 
c s = in Eq.2. Additionally, the entropy at the transi- 
tion temperature determined from the specific heat in the 
model is within a few percent of that of the normal fluid 
at T c , as is required for a second order thermodynamic 
transition. Although Eq.2 is highly phenomcnological, 
nonetheless it might be a useful guide over a wider range 
of temperature than we have explored. It gives a low 
temperature limit for the density of states, relative to the 
normal fluid, to be simply proportional to the pressure 
dependent coherence length, a£o(P)A/Vi. It would be 
interesting to extend heat capacity experiments to lower 
temperatures for a direct measurement of the density of 
states of surface bound states. 

In conclusion, we have used a high resolution method 
to determine the heat capacity of 3 He in the presence of 
a silver surface. We distinguish two different contribu- 
tions; one from the bulk superfluid phase and the other 
from the helium near the silver surface. We have con- 
structed a model based on low energy contributions to the 
density of states associated with Andreev bound states 
of 3 He quasiparticles that scatter from the surface. We 
have found that the surface heat capacity has a temper- 
ature and pressure dependence given by the normal fluid 
specific heat and the bulk 3 He coherence length. Fur- 
ther, we determine that the spatial extent of the bound 
state region is one half of the bulk 3 He coherence length. 
Our confirmation of the existence of surface bound states 
from measurement of the heat capacity supports recent 
results from surface sensitive measurements of the trans- 
verse acoustic impedance by Aoki et raZ.[lfi| and theoret- 
ical calculations of superfluid 3 He in slabs by Vorontsov 
and Sauls 0j- 
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